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We investigate the effect of different edge types on the statistical properties of both the energy 
spectrum of closed graphene billiards and the conductance of open graphene cavities in the semi- 
classical limit. To this end, we use the semiclassical Green's function for ballistic graphene flakes 
that we have derived in Reference Q]. First we study the spectral two point correlation function, or 
more precisely its Fourier transform the spectral form factor, starting from the graphene version 
of Gutzwiller's trace formula for the oscillating part of the density of states. We calculate the two 
leading order contributions to the spectral form factor, paying particular attention to the influence 
of the edge characteristics of the system. Then we consider transport properties of open graphene 
cavities. We derive generic analytical expressions for the classical conductance, the weak localization 
correction, the size of the universal conductance fluctuations and the shot noise power of a ballistic 
graphene cavity. Again we focus on the effects of the edge structure. For both, the conductance 
and the spectral form factor, we find that edge induced pseudospin interference affects the results 
significantly. In particular intervalley coupling mediated through scattering from armchair edges 
is the key mechanism that governs the coherent quantum interference effects in ballistic graphene 
cavities. 

PACS numbers: 72.80.Vp, 73.22.Pr, 05.45.Mt, 73.23.Ad, 03.65.Sq 



I. INTRODUCTION 

Nearly twenty- five years ago, Berry and Mondragon^ 
introduced the 'neutrino billiard', a Dirac Hamiltonian 
describing a massless spin one-half particle in a plane 
with lateral confinement, as a fictitious, conceptually 
simple quantum system which exhibits time-reversal 
symmetry (TRS) breaking without magnetic fields. They 
demonstrated that the energy eigenvalues of such a Dirac 
billiard are statistically distributed according to the cor- 
responding Gaussian unitary ensemble (GUE) of random 
matrix theory (RMT), i.e. the ensemble of random her- 
mitian matrices without TRS, if the quantum system 
possesses a classical counterpart that is chaotic. 

Ballistic cavities built from monolayer graphene (for re- 
cent reviews see Refs.^-d) with Fermi energy close to the 
Dirac point, have been proposed as realizations of such a 
neutrino billiard^—. In fact due to the coexistence of the 
Dirac points associated with the two independent valleys 
in the graphene band structure, graphene-based billiards 
represent two copies of a Dirac billiard mutually coupled 
through intervalley scattering. For vanishing intervalley 
coupling, the entire graphene Hamiltonian is composed 
of two decoupled sectors with degenerate sets of eigenval- 
ues each set obeying GUE statistics just like the neutrino 
billiard. However, for sufficiently strong intervalley cou- 
pling, one expects an eigenvalue statistics according to 
the Gaussian orthogonal ensemble (GOE), because the 
entire graphene Hamiltonian is time-reversal invariant, 
i. e. it is a real symmetric matrix in the absence of ad- 
ditional magnetic fields. Therefore intervalley coupling 



is an effective symmetry restoring parameter that allows 
tuning from GUE to GOE behavior in graphene. 

In disordered graphene samples, short ranged scatter- 
ing potentials or vacancies cause intervalley scattering. 
Depending on the corresponding intervalley scattering 
time and other relevant inverse scattering rates of the 
system (c. g. the phase coherence time), the magnetocon- 
ductance of disordered bulk graphene is expected to show 
signatures of weak localization or weak antilocalizatioiA 
This theoretical picture has been confirmed experi- 
mentally to some extent by the measurement of weak 
localization^— and weak antilocalization-^ in disor- 
dered graphene. 

In contrast, in ballistic graphene nanostructures (with 
clastic mean free paths on scales of several microns^) , it 
is the edges that are responsible for intervalley scatter- 
ing. However, not all types of edges arc effective inter- 
valley scattcrcrs: While scattering of quasiparticlcs from 
a zigzag or an infinite mass type edge is an intravalley 
scattering process, armchair edges couple quasiparticle 
states from both Dirac points. Recent experiments on 
nanostructured graphene^ have revealed that the mag- 
netoconductance in such samples cannot be consistently 
described using the bulk theory of Ref.|^ and signatures of 
ballistic transport were found to be non-negligible. While 
several numerical studies have investigated quantum in- 
terference effects in spectral and transport properties of 
ballistic graphene nanostructures^—, so far there is no 
analytical theory that is capable of treating quantum in- 
terference effects in arbitrarily shaped graphene nanos- 
tructures. In this work, we provide a semiclassical the- 
ory for the conductance of ballistic graphene structures 
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as well as for spectral correlations, particularly focusing 
on the effects of different edge types. Our semiclassical 
approach requires L> A E with A E the Fermi wavelength 
and L the system size, while at the same time the ener- 
gies should be compatible with the assumption of a linear 
dispersion law in the effective Dirac theory for graphene, 
i. c. X E <C a, with the graphene lattice constant a. Fortu- 
nately, both conditions are well fulfilled for most meso- 
scopic graphene structures: For example, for a system 
size of order L ~ 100 ran, there are thousands of levels 
in the linear dispersion regime, where the Dirac equation 
is valid. 

This is the second paper of a two-paper series on edge 
effects in graphene nanostructures. In the first paper—, 
referred to as paper I in the following, we derived an ex- 
pression for the single-particle Green's function of a bal- 
listic or weakly disordered graphene structure in terms 
of multiple reflections from the system boundary. Using 
this expansion, we could account for different types of 
boundary conditions (zigzag, armchair, infinite mass or 
combinations of those) , while the propagation inside the 
system was treated on the level of the effective Dirac the- 
ory. Furthermore, we obtained expressions for the mean 
density of states (DOS) valid in the semiclassical limit, 
highlighting the role of edge states at zigzag boundary 
segments, as well as for the oscillating part of the DOS in 
terms of periodic orbits for classically chaotic and certain 
intcgrable graphene cavities. Graphene edge phenomena 
and the effects from sublattice and valley pseudospin dy- 
namics enter into those trace formulae through (traces 
over) pseudospin propagators evolving along the classical 
orbits that otherwise are the same as in a corresponding 
non-relativistic billiard. 

Here we consider the spectral statistics of closed 
chaotic graphene cavities as well as the conductance of 
open graphene cavities coupled to (two) leads. We ad- 
dress the question as to how their universal properties 
are governed by the underlying Dirac-type dynamics to- 
gether with the graphene-specific edge effects. In particu- 
lar, we study signatures of the afore mentioned crossover, 
mediated through intervalley scattering, between unitary 
and orthogonal symmetry classes in spectral and trans- 
port observablcs. We are thereby able to give semiclassi- 
cal explanations and analytical expressions for results ob- 
tained earlier in Ref.[l9| mainly numerically and in terms 
of symmetry arguments. There it was shown that inter- 
valley coupling (mediated through scattering at armchair 
edges) acts differently on spectral and transport prop- 
erties of classically chaotic graphene cavities: While in 
quantum transport a crossover from unitary to orthogo- 
nal behavior should be observable with increasing inter- 
valley coupling, spectral statistics probing the scale of the 
mean level spacing is predominantly governed by the or- 
thogonal symmetry class, even for rather weak armchair 
scattering. Here we study the spectral statistics at en- 
ergy scales larger than the mean level spacing, where the 
unitary symmetry should be observable even at moderate 
intervalley scattering. 



To this end we start from the semiclassical expressions 
for the Green's function and DOS in terms of interfer- 
ing classical trajectories, derived in paper I, which we 
introduce in Sec.|TTJ In Sec.|HT] we address spectral fluc- 
tuations, more specifically we derive semiclassical expres- 
sions for the spectral form factor, the Fourier transform of 
the spectral two-point correlator characterizing spectral 
statistics. We find that the leading order and the next 
to leading order contributions are strongly influenced by 
the edges. The total amount of armchair edges is the rel- 
evant quantity that determines the size of graphene spe- 
cific correlations: Our results suggest that evidence for 
partial (effective) breaking of the time reversal symme- 
try (TRS) is visible in the spectral correlations at scales 
of many mean level spacings. The reader with main in- 
terest in transport can also skip Sec. lIIII and directly go 
to Sec. lIVi which is devoted to mesoscopic phenomena 
in coherent transport. We study imprints of the edges 
and the chaotic dynamics in confined graphene systems 
on weak localization, conductance fluctuations, and shot 
noise. 

We start from an exact expression for the Green's 
function in terms of multiple scatterings from the 
boundaryi22r— and generalize advanced semiclassical 
techniques, introduced in Refs.[23| and [24| for graphene, 
accounting for classical trajectory correlations, which 
have proven essential for understanding spectral statistics 
and quantum phenomena in ballistic transport through 
chaotic Schrodinger billiards. In this way we derive an- 
alytical expressions for the observablcs mentioned above 
that allow us to predict the edge dependence of several 
quantities in the semiclassical regime. While the details 
of the edge structure are crucial for the form factor, the 
weak localization correction and the conductance fluc- 
tuations, it turns out that shot noise is not affected to 
leading order in the inverse channel numbers of the leads. 
We finally conclude and give an outlook in Sec.fV] and 
collect longer derivations in three appendices. 



II. GREEN'S FUNCTION AND ANTIUNITARY 
SYMMETRIES 

A. Semiclassical Green's function for graphene 
billiards 

In Sees. Mil and HVI wc employ the Green's function in 
order to study correlations in the spectra and the con- 
ductance of graphene flakes in the semiclassical limit. In 
paper I we derived an exact expression for the Green's 
function of a graphene flake in terms of a multiple reflec- 
tion expansion, where each term in this expansion cor- 
responds to the sum of Feynman paths with a specific 
number of reflections at the system boundary. In this 
section we repeat our main results that will be useful for 
the present paper. 

For a confined ballistic graphene structure quasiparti- 
cle dynamics is described by the effective Dirac Hamilto- 
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with appropriate boundary conditions at the edges. Here 
the {tJj} denote Pauli matrices in sublattice pscudospin 
space and Pauli matrices in valley-spin space are repre- 
sented by {t^}, while oo and To are unit matrices acting 
on the corresponding spin space. Then the full Green's 
function is given in the semiclassical limit by (cf. Eq. (78) 
in I) 




G(x,x') 



( JSy/h+ifj-yTT/2 



(2) 
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Here the sum runs over all classical orbits 7 leading from 
point x' to x. These classical orbits consist of straight 
lines that are connected by specular reflections at the 
boundary, so that the condition of least action is satis- 
fied. This situation is similar to billiards with spin orbit 
interaction (SOI) where the entire effect of the SOI is con- 
tained in the spin matrix with the same trajectories^—. 
In Eq. @ Sj = hk E L 1 and // 7 are the classical action 
(k E and L 1 are the Fermi momentum and the length, 
respectively) and the number of conjugate points for the 
classical orbit 7, and 



D 7 = 



1 



dx± 



9p'± 



-1/2 



(3) 



denotes an element of the stability matrix of the path 7. 
Here p' ± and x± are the components of the initial momen- 
tum and final position perpendicular to the trajectory. In 
Eq. ©, 



(4) 



is the pseudospin propagator along the orbit 7, with 
u a x > = ct ■ (a 1 — x'). The product in Eq. runs 
over all iV 7 reflections that occur at boundary points a.i 
along the orbit 7, with the edge dependent pseudospin 
rotations for reflections at zigzag (zz), armchair (ac), and 
infinite mass (im) type edges 



Ki = ±i 



^ „i2Kxn 



2> e i 



for zz , 
for ac , 
for im . 



Here Oi denotes the reflection angle at the z-th bounce 
with the boundary, while ti and Xi are the direction of 
the tangent to the boundary and the x-coordinate of the 
reflection point cti, respectively (see Fig.[T]). Further we 
define for a given vector v the Pauli matrix <j v = <j ■ v. 
K = 4ir/3a is the distance of the Dirac points from the 
T-point of the Brillouin zone. For reflections at infinite 
mass type edges, the sign in Eq. ([5]) is equal to the sign 
of the mass potential outside of the system^. For zigzag 
edges the sign is determined by the type of sublattice at 



FIG. 1. Scheme of two classical reflections along an orbit 7. 
The normal directions at the reflection points on and au+i 
are denoted by unit vectors rii and Tii+i, respectively. The 
tangential directions are ti = rii x z and ti+i 
while the reflection angles are Qi and respectively. 



the zigzag edge. For an A-edge, the upper sign is valid 
and for a B-edge the lower sign. For armchair edges, the 
upper sign is valid when the order of the atoms within 
each dimer is A-B along the direction of ti, and the lower 
sign is valid for B-A ordering. 

Equations ([2] - [5]) specify the contribution of a given 
classical orbit 7 to the semiclassical Green's function. 
Note that the dynamical part is identical to the semi- 
classical Green's function of a Schrodinger system, while 
it is the pseudospin propagator K-y, through which the 
graphene specific physics enters. This relation allows 
us to use many results from the semiclassical theory of 
Schrodinger billiards. 



B. Antiunitary symmetries of the effective theory 

In order to compare the results from our semiclassical 
theory with universal RMT predictions, we need to con- 
sider the relevant antiunitary symmetries of the Hamil- 
tonian H, Eq. (fTJ). An antiunitary operator is given by 
the complex conjugation followed by a unitary operator. 
We assume that all spatial symmetries, such as invari- 
ance under reflection or rotation, are broken. Thus the 
only relevant unitary operations are rotations in the pseu- 
dospin spaces. As mentioned above we adapt here the 
valley-isotropic basis, in which the bulk Hamiltonian is 



(5) invariant under the transformations 7^ HTi with=* 



•nC, ie{0,x,y}. 



(6) 



1 while T? 



Furthermore we note that Ty — ± wmiu ; 

the overall TRS that connects states from 



T 2 = -1. 

That means T y , 
opposite valleys, drives the system into the orthogonal 
symmetry class. This symmetry is not broken by any of 
the boundaries we consider. However, we will see that it 
is only relevant if intervalley scattering is present. In the 
absence of intervalley scattering, the effective intravalley 
TRS 7o gives rise to the symplcctic universality class. 
However, for the boundary conditions that do not couple 
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the valleys (infinite mass, zigzag), 7o is not preserved, 
and thus it is irrelevant for ballistic systems. Finally T x 
is also a symplectic symmetry. It is broken by intcrvalley 
scattering armchair edges and left intact by zigzag or 
infinite mass type edges. For preserved valleys it thus 
ensures the (Kramers) degeneracy of the two valleys. 

We note that, since the bulk Hamiltonian has a 
particle-hole symmetry (<j x H*o~ x = —H), the chiral uni- 
versality classes are relevant at zero energy. In this pa- 
per, however, we consider Fermi energies away from zero, 
where the normal and the chiral universality classes lead 
to the same results 2 ^. 

Sections |TlT]on spectral statistics and Sec.|lV]on trans- 
port can both be read independently. We begin with 
the analysis of spectral statistics which is conceptually 
slightly simpler. 

III. SPECTRAL FORM FACTOR OF CHAOTIC 
GRAPHENE BILLIARDS 

A. Spectral form factor: 
Definition and random matrix results 

Quite generally, the DOS can be split into a smooth 
part, i. e. the mean DOS p, and an oscillating part p osc , 

p(k E ) = p{k E ) + p osc (k E ) . (7) 

While p contains information about the coarse structure 
of the spectrum, details about the level distribution are 
in p osc . The mean DOS p and the trace formulae for the 
oscillating DOS p sc were addressed in I. Here we use the 
trace formula for p osc to study statistical properties of 
the spectral fluctuations. To this end, we focus on the 
spectral two-point correlator 

R(v) = ^7T-^(Posc(fc E + ?7/2)posc(fc E - v/ 2 ))k E , (8) 
P{k E ) 

where ( }f- B denotes averaging over a (classically) small 
window of the Fermi energy hv F k E or, respectively, on 
its Fourier transform the spectral form factor 

oo 

F(t) = J driFt^y 2 ^^ . (9) 

— oo 

Here t denotes the time in units of the Heisenberg time 
T H = 2irp(k E )/v F . In the limit t -> 0, RMT predicts the 
universal expression^ 

( 2t-2t 2 + (D(t 3 ) forGOE, 
F(t) = \ t for GUE , (10) 

[ \t+ \t 2 +0{t 3 ) for GSE, 

which is expected to be valid for chaotic systems. Sys- 
tems with a TRS T and T 2 = 1 (e.g. a free electron 
without magnetic field) belong to the orthogonal sym- 
metry class and systems without TRS (e.g. an electron 



in a magnetic field) to the unitary class. Systems with 
a TRS T and T 2 = —1 (e.g. a spin 1/2 particle with 
broken spin rotational symmetry) are members of the 
symplectic class. However, as mentioned above, for bal- 
listic graphene flakes the edges break the symmetry 7o 
that would lead to the symplectic universality class. In 
the presence of spin-orbit interaction, which we neglect 
in this paper, the symplectic class could be relevant. 



B. The spectral form factor: 
Semiclassical evaluation 

1. Gutzwiller trace formula for graphene 

For chaotic Schrodingcr billiards, it is possibl o 23 ' 31 i 32 
to connect the properties of classical orbits to F(t). Using 
this relation the universal RMT form factor (fTOj) has been 
derived semiclassically using Gutzwiller's trace formula^ 3 - 
for Schrbdinger systems. In the following we will extend 
this semiclassical theory to describe the spectral correla- 
tions of chaotic graphene systems. 

We start by the following extension of Gutzwiller's 
trace formula (see Eq. (103) in I) to graphene: 

Posc(k E ) = ^Re V A 7 Tr(A' 7 )e^ L - . (11) 

7 

Here the sum runs over infinitely many periodic classical 
orbits 7. The classical amplitudes A~ i depend on the 
period, the stability and the number of conjugated points 
of the corresponding orbit 3 - 3 -. As mentioned above, the 
trace over the pseudospin propagator K 7 accounts for 
the graphene features. 

From Eqs. ([HE]) we know that the valley part of K 7 
contains an off-diagonal factor 

TKi — Ty e (12) 

for every reflection at an armchair edge. This implies that 
the trace vanishes, Tr K 7 = 0, unless the total number 
of ac reflections iV ac along 7 is even. If iV ac is even, 
the pseudospin factor we have for an orbit with a total 
number of N = N 7 reflections is (see App. lA lj ) 

Tr K 7 = 4/ 7 cos (0 7 + Ntt/2) 

x cos(2is:A 7 + ?9 7 + Ntt/2) , (13) 

where / 7 is a phase factor that depends on the exact 
sequence of ac, zz, and im type reflections, 

N 

7 = $> (14) 

i=l 

is the sum over all reflection angles 0i along the orbit 7, 

JVac/2 

A 7 = ^2 (x 2 i-i - X 2 i) (15) 

i=l 
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FIG. 2. Scheme of the orbit pairs that contribute to the spec- 
tral form factor, Eq. (|17[1 , to leading order beyond the diag- 
onal approximation. The orbit 7 X intersects itself with an 
angle e, while the orbit 7 avoids this crossing. Apart from 
this encounter region, the orbits follow each other exponen- 
tially closely, with the same propagation direction in the left 
loop L and opposite direction in the right loop R. Note that 
the real trajectories in a billiard consist of a series of many 
straight pieces and classical reflections at the system bound- 
ary between them, as shown in Fig.[T] 

measures the differences in the ^-coordinate of pairs of 
subsequent ac reflection points, and 

^ = Y,(-l) Si #i- (16) 
i=i 

In Eq. ([To 7 ]) the sum is restricted to the N zz reflections at 
zz edges along 7, and s, is the number of ac reflections 
that occur after the specific zz reflection i. Furthermore 
$i = +9{ for reflections at ^4-edges and = — Oi for 
reflections at S-edges. We note that TrA' 7 = Tr K~-i 
where 7 -1 and 7 are time reversed partners in the clas- 
sical sense, i. e. they are identical up to the direction of 
movement. 



2. Products over phase- carrying paths 

We now evaluate the spectral form factor for graphenc 
semiclassically. We insert Eq. (|TT|) into the definition of 
F(t) and obtain the approximate expression 

(17) 

where T 7 = L 7 /u F is the period of an orbit, T = £Th, 
and the short notation for the pseudospin traces 

Zy.y =TrKr i TcKl l , . (18) 

The action difference between the orbits is given by 
SS'y^' —— tik E (^L^ _Zv-y' ) . 

If (7,7') is an arbitrary, uncorrclated pair of orbits, 
e *ss y ^,/h j g g ener i ca Uy a rapidly oscillating function of 
k E in the scmiclassical limit. Thus after k E averaging, 
the corresponding contribution to F vanishes and the 
relevant contributions to F come from classes of orbit 
pairs with small or even vanishing action difference. That 
means the actions of the orbits have to be classically cor- 
related. The most obvious class of correlated orbit pairs 
are captured by the so-called diagonal approximation, 



where only equal action pairs 7' = 7 (and 7' = 7 _1 ) 
are considered. These orbit pairs give rise to the leading 
order (~ t) terms in the expansion ([T0| for F(t)Q. The 
second order terms (~ t 2 ) in the GOE case are connected 
to another class of periodic orbit pairs: loop contribu- 
tions that are sketched in Fig. [2 These two orbits follow 
each other closely for most of the time but one of them 
has a self-crossing (7 X ) while the other one 'avoids' this 
crossing (7). Therefore the propagation directions are 
the same in one part of the paths (L) and opposite in the 
other part (R)2S. 

In the following we study the diagonal contribution 
and the off-diagonal loop corrections (Fig. [2} to F(t) for 
graphenc billiards. We find that both contributions de- 
pend on the structure of the edges due to the interplay 
between the boundary conditions and the (effective) TRS 
of graphenc flakes. 

3. Diagonal contribution 

For the diagonal terms in the double sum (jTT|) . the 
spectral form factor reads (including pairs 7' = 7 and 
7' = 7-1) 

7 

Equation (jT5)l yields for the product of traces (|T51) 

Z-y == Z-y^ == Zry™ — 1 (^0) 

= 16cos 2 (6' 7 + iV7r/2)cos 2 (2 J fCA 7 + 79 7 + iV7r/2) , 

provided that the total number of ac reflections N ac along 
7 is even, otherwise Z 7 — 0. Since 7 is a periodic orbit, 
the total rotation angle is 2ttw, where w is the integer 
winding number. Therefore we can write 

7 = (7V+-iV_)|- W 7r, (21) 

where N + and iV_ count the reflections with positive 
or negative reflection angles, respectively. This leads to 
cos 2 (6> 7 + Nir/2) = 1 and thus for even A^ ac 

Z 7 = 16 cos 2 (2XA 7 + t? 7 + Nn/2) . (22) 

a. No ac scattering - valley conservation 

For pedagogical reasons, we first consider the simpler 
case of billiards whose boundaries consist only of a combi- 
nation of zz and im type edges. In other words we assume 
that iVac = for all orbits, so that there is no interval- 
ley coupling at all. In order to perform the orbit sum in 
Eq. (fTT|) we need to average Z 1 accordingly. To perform 
this averaging, we first note that classically chaotic bil- 
liards exhibit ergodic dynamics, thus a long orbit will hit 
all boundary points with the same probability. Therefore 
we average Z 7 over the possible edge types or boundary 
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conditions, respectively. Since we assume that no ac edge 
are present A 7 = in Eq. (f2"2"j) . If no zz edges are present 
cither, also $ 7 = 0, otherwise $ 7 is randomly distributed 
between zero and 2tt. In both cases the squared cosine 
in (|22p is on average 1/2, since for long orbits even and 
odd N are equally probable. That means we obtain 



(23) 



Now we can pull (Z 7 ) out of the sum in (fH)]) 
and are left with the evaluation of the form fac- 
tor as in the Schrodinger case. By employing the 
Hannay-Ozorio dc Almeida sum rule^i we obtain^ 

F D (t)^^-D(T)\A(T)\ 2 = 4t. (24) 

According to Eq. (fTU|) this result is four times the GUE 
prediction for the linear term of F(t). To understand 
the factor of four we recall the antiunitary symmetries 
of the problem. The physical TRS T y would lead to 
the orthogonal symmetry class. However, since there 
is no intervalley scattering and 7~ y connects states from 
different valleys, it is basically irrelevant. The effective 
intravalley TRS To, which would lead to a symplectic 
ensemble, is broken by both zigzag and infinite mass 
type edges. Hence we end up with two unitary subsys- 
tems that are Kramers' degenerate partners due to the 
remaining symmetry T x . Equations (|5J) and © yield 
F{t) = 2.Fi(2t) in this case, where F\ is the form factor 
of the non-degenerate subsystem. The loop contribution 
(see below) agrees with this picture. 

b. ac billiards - complete valley mixing 

We consider the opposite limit and assume that all 
reflections happen at ac edges. The term KA 1 can only 
take discrete values 



if A 7 = —n n £ Z . 
3 



(25) 



For a generic boundary, n is random and thus the aver- 
age of the pseudospin contribution (|22|) reads (Z 7 ) = 4. 
Consequently we obtain the diagonal contribution to the 
spectral form factor 



F D (t)=2t, 



(26) 



in agreement with the RMT prediction for the or- 
thogonal symmetry class. The orthogonal symmetry 
is expected because the intervalley mixing due to the 
armchair scattering restores the TRS described by %^ ■ 

c. Mixed boundaries - the generic case 

In the realistic situation of mixed boundaries, 
i. e. boundaries that consist of both, valley conserving 
edges and valley mixing edges, we obtain, based on our 
earlier discussion, that 



where P° c (T) denotes the average probability that an or- 
bit with period T hits armchair edges an even number of 
times (including zero times). Assuming chaotic dynam- 
ics and that the total length of ac edges of the billiard 
W ac = |<9V ac | is small compared to the total boundary 
length, the ac scattering rate is approximately given by 
(see Ref-Hil and references therein) 



1 

t7 c 



ttA 



(28) 



The armchair scattering time T ac is the typical time be- 
tween two reflections at ac type edges. In this case P ac 
is approximately given by 



P' 



^H( 



1 + e 



-2T/T a , 



(29) 



With that we can evaluate the orbit sum ((IT)) as before 
to obtain 



F D (t) = 4tP a c c (T) = 2t(l + e - 2tTH/T > 



(30) 



This formula describes the crossover between the two lim- 
iting cases of completely decoupled and mixed valleys, 
respectively. For the asymptotic limits Th/T^c — > and 
Tn/T ac -^oowe recover Eqs. (|24)l and ([26]) . respectively. 

An advantage of Eq. (|3"0|) is that the crossover param- 
eter does not have to be introduced phcnomcnologically, 
but can be directly related to the microscopic proper- 
ties of the system, namely the total amount of (interval- 
ley scattering) ac edges. The ratio of the relevant time 
scales is directly proportional to the total amount of ac 
boundary pieces, namely: 



T 

± ac 



2fc E W a , 



(31) 



(27) 



One can understand this with the following construction: 
if the intervalley relaxation is modeled in the spirit of 
Biittiker as intervalley relaxing lead with width W ac , then 
Tn/T ac is proportional to the number channels of the 
intervalley relaxing lead. 

Figure [3]a) shows F D (t)/t for small t = T/T H and the 
parametric crossover between the limiting RMT sym- 
metry classes, namely GOE [lowest line, Eq. ([26]) ] and 
4xGUE [uppermost line, Eq. (|2"4"|)]. In the general case 
of finite Tn/T ac , Fp is not purely linear as in both lim- 
iting cases. Note that even for a considerable amount 
of ac scattering (consider e. g. the (turquoise) curve for 
Tu/T ac — 10) we still find rather strong deviations from 
the GOE statistics. Panel b) shows the crossover of F(t) 
as a function of Tn/T ac . We point out that this crossover 
from 4xGUE to GOE is rather particular: Usually tran- 
sitions between GUE and GOE result from a symmetry 
breaking mechanism, e. g. due to a magnetic field. Here 
the armchair edges and the resulting intervalley scatter- 
ing act as a symmetry restoring mechanism. Moreover 
the correlations in the unitary limit are weaker than in 
the orthogonal case, in contrast to the case of the usual 
GOE-GUE transition. 
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a) 




0.02 0.04 0.06 0.08 0.10 



b) 




10 20 30 40 50 

T H /T ac 

FIG. 3. a) Rescaled diagonal contribution Fn(t)/t to the 
spectral form factor as a function of t = T/Th for a generic 
graphene cavity for different values of the intervalley scat- 
tering strength. From top to bottom: Tu/T^c = (black), 
1 (red), 3 (blue), 10 (turquoise), 50 (orange) and oo (green). 
While in the limiting cases of TH/T ac — ► and Tn/T ac — s> oo 
the i-dependence of Fo is purely linear, this is not generally 
the case. For all finite values of TH/Tac, Fd shows signatures 
of unitary correlations at small t. b) From bottom to top, 
F D as a function of T H /T ac at t = 0.02 (black), 0.04 (red), 
0.06 (blue), 0.08 (turquoise) and 0.1 (orange). It is an expo- 
nential crossover from 4 times the GUE to the GOE result. 



In a mesoscopic graphene quantum dot the crossover 
parameter Th/T^c can be tuned by changing the energy. 
Since our theory is valid in a wide range of k E , the tran- 
sition in Fig.[3]b) should be accessible in a real system. 
Consider for example a structure with a typical length 
scale of 1000 a, then we expect our theory to be valid for 
0.01 < k E a < 0.2, approximately. Assuming W ac « 250 a 
this leads to 5 < T H /T ac < 100. 



4- Off-diagonal (loop) contributions 

Now we proceed with the calculation of the correction 
terms (order t 2 ) to the form factor. At this order it is 
the loop pairs (7, 7x) depicted in Fig. [51 that contribute 
to F(t). We denote these contributions by Fl- Along 
the two loops L and R the two trajectories are exponen- 



tially close to each other. For vanishing Ehrenfcst time, 
we assume that the sequences of reflections along 7 and 
7x are exactly equal in the L-part while they are equal 
but with opposite order in the i?-part. Orbit pairs that 
differ in their number of reflections are of measure zero 
in the semiclassical limit. We start by calculating the 
pscudospin propagators. All quantities will be labeled 
correspondingly for the loops L, R. For the total reflec- 
tion angles we have 

d-y = 0L + 0R j #7 X = &L — 0R • (32) 

As we show in App. lA 2l the expression (flU)) for the traces 
leads to 

Z 7i7x = (-l) w " +1 16cos(2ifA 7 + ?? 7 + A%/2) 

x cos (2KA 7x + + Ntt/2) (33) 
if ./Vac is even for both orbits, and Z ltly = otherwise. 

a. No ac scattering - valley conservation 

When there is no ac scattering, the valleys are uncou- 
pled and we have 

Z 77x = (-l) Ar « +1 16cos(tf 7 + Ntt/2) cos (tf 7x + Ntt/2) . 

(34) 

We decompose the angles # into $l ± $r and thus obtain 

cos (i? 7 + Ntt/2) cos (t? 7x + Ntt/2) (35) 
= cos 2 + Ntt/2) cos 2 (& r ) 
- sin 2 {§ L + Ntt/2) sin 2 {d R ) . 

The averaged pseudospin contribution to Fl vanishes 

(Z 7 , 7x )=0, (36) 

because the angles •Ql and have to be averaged sep- 
arately, and thus the cosine terms and the sine terms 
cancel on average. This result holds also if zz edges are 
absent and all reflections happen at im edges. In the 
latter case t? 7 = ^ 7x = and therefore ^ 7l7x is either 
identically zero (odd N) or Z 7 . 7x ~ (—1)^ (even N), 
leading again to Eq. (|3"6")l . because for long orbits there 
arc equally many orbits with even and odd Nr£. That 
means there is no loop contribution to the spectral form 
factor, 

F L (t) = 0, (37) 
if there are no intervalley scattering ac edges. 

b. ac billiards - complete valley mixing 

In the opposite limit of dominant ac scattering, 
Eq. (33]) gives for even N = N ac 

Z 7i7x = (-l) JVK+1 16cos(2is:A 7 )cos(2 J fs:A 7x ) . (38) 

To average this, we have to distinguish between two dif- 
ferent cases. First we assume that N R , the number of 
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reflections in the i?-part of 7, is even. Then Al and A R 
are well defined and we have 



A 7 = Al + A R , A 7 



Ar - A, 



(39) 



From the trigonometric relation analog to Eq. (|35|) it fol- 
lows that these orbits do not contribute on average. On 
the other hand if N R is odd, it is straight-forward to show 
that 



cos(2A'A 7x ) = cos(2A"A 7 ) 
and therefore for even N 



J 7,7x 



16cos 2 (2A"A 7 ). 



(41) 



Since only orbits with even TV and odd N R contribute we 
obtain for the average 



Z 7 , 7x ) = 4 (cos 2 (2 K A 7 )) = 2, 



(42) 



that is the graphene specific pseudospin trace enters 
merely with an overall prefactor 2 to the spectral form 
factor. The rest of the calculation is identical to eval- 
uation of the semiclassical form factor in Ref.[23j. The 
full contribution of pairs of orbits (7, 7x) is obtained by 
summing over all 7 and for each 7 counting the number 
of self encounters, i. e. self intersections with a (small) 
crossing angle e. We introduce P(e, T 7 ) as a measure for 
the density of self crossings of an orbit of length Ty 2 ^, 



T-T min ( £ ) 



P(e,T)«^ J dT R [T-T R }. 



(43) 



Here 



T min (e) 



T m i n (e) = --ln(ce) 



(44) 



with A the average Lyapunov exponent and c a constant, 
is the minimal time to form a closed loop. Then the an- 
gular integral over the exponents of the phase differences 
5S 7 . lx /h is computed, leading to— 



71 

Re J dee lSS '<--<* /h P(e,T 1 ) = -t 



(45) 



as for a Schrodingcr billiard^. Evaluating the remaining 
orbit sum similar to Eq. (|19[) yields 



(46) 



F L {t) = -2t\ 



in agreement with the RMT prediction [cf. Eq. (|10p] for 
the GOE. 



c. Mixed boundaries - the generic case 

We now consider the general case, where all three types 
of edges - ac, zz, and im - are present. From the previ- 
ous discussion we know that the contribution to Fl is 
exclusively due to such classical orbits 7 that undergo 
odd numbers of ac reflections along both, the L (TV aCj £) 
and the R (TV aCl i{) parts of 7. For these orbits we ob- 
tain (Z 7 7x ) = 8 also in the presence of zz and im type 
reflections, while (Z lt7x ) = for the remaining orbits. 



(40) Therefore we can write, analogous to Eq. ([45 



F L (t) = 8iRe / dee 



iSS-, 



P(e,T 7 ) (47) 



with P(e, T 7 ) the corresponding density of self crossings. 
In order to determine the density of self crossings of orbits 
that fulfill the condition of odd N aCtL / R , we start from 
the unrestricted density of loop pairs (|4"3"]l . To account 
for the fact that only the subset of loop pairs with odd 



N., 



,L/R 



contributes, we have to modify P to obtain 



vie 



T-T m 



P(e,T)«^ / dT R [T~T R ] (48) 

T mln (e) 

xP° c {T R - T min ( £ )]P°[T - T R - T min (e)] 

with the approximate probability to hit ac edges an odd 
number of times during a time T 



K C [T) = \ ( 



-2T/T a , 



(49) 



Using Eq. (|49[) when evaluating the integral in Eq. 
yields 



P( £ ,T 7 ) 



vjeT 
wA 8 

-2T min (e) 



T - Tac + (T + T ac )e~ 2T / T - (50) 
2T~ 



1-1 



D -2T/T ac 



Here we assumed that (in the RMT limit) T m i n (e) is much 
shorter than the other time scales T and T ac and there- 
fore neglected terms of higher order in T m ; n (e). The e- 
independent terms of P do not contribute 2 ^ to the real 
part of the integral in Eq. (|4T|) . so that we get 



F L {t) 



-2e 



[ 1 + 2 HA e -2 t T H /T a< 



(51) 



We recover the asymptotic limits (f3"T)) and (|4^|) for 
Tn/T ac -> and T H /T a c -> 00, respectively. 

Figure |4]a) shows —Fhifyjt 2 for various effective ac 
scattering strengths Tn/T ac in the range of small t, where 
Fl is the dominant off-diagonal contribution to the form 
factor. The parameter T^/T ac controls the crossover be- 
tween the two RMT limits: F L = for T H /T ac -> and 
F L = -2t 2 for Tft/T ac — >• 00. In panel b) we plot the 
crossover of —Fl(£) as a function of T^/T ac . 
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a) 




10 20 30 40 50 

T H /T ac 

FIG. 4. a) Absolute value of the rescaled leading off-diagonal 
contribution, —FL(t)/t 2 , to the form factor as a function of 
t = T/Th for different values (from bottom to top) of the 
ac scattering strength Th/T&c = (black), 5 (red), 10 (blue), 
20 (turquoise), 50 (orange) and oo (green). At Tn/T ac — Fl 
is vanishing, in agreement with GUE, while for larger values 
the curve approaches the quadratic behavior of the GOE pre- 
diction, b) Exponential crossover of — Fl as a function of 
Tk/Tkc (from bottom to top) at t = 0.02 (black), 0.04 (red), 
0.06 (blue), 0.08 (turquoise) and 0.1 (orange). 



Equations (f3"U)) and (|5Tj) are the main results of this sec- 
tion. They describe the spectral correlations of chaotic 
graphene billiards with a boundary that consists of an ar- 
bitrary combination of ac, zz and im type edges. The to- 
tal amount of intervalley scattering at armchair edges sets 
the time scale T ac that represents the control parameter 
for a 4xGUE to GOE crossover. In App.lBlwc generalize 
our results by incorporating direct TRS breaking, e. g. by 
including a small magnetic flux $ through the billiard. 
Due to the flux induced Aharonov-Bohm type phase dif- 
ferences, this causes a breaking of both, the TRS T y and 
the valley symmetry T x leading to another crossover. In 
the asymptotic cases of vanishing ac edges and complete 
valley mixing, we find transitions 4xGUE— s- 2xGUE and 
GOE-s- GUE, respectively. 



IV. TRANSPORT THROUGH OPEN 
GRAPHENE CAVITIES 

The crossover between the different effective universal 
symmetry classes discussed for spectral statistics in the 
proceeding chapter is also reflected in quantum transport 
properties that are experimentally more directly accessi- 
ble. In this part of the paper we focus on the two termi- 
nal conductance of open graphene cavities. Our starting 
point is the linear response expression for the conduc- 
tance from lead b to lead a in terms of the Green's func- 
tion of the cavity^ 

g a b = - j dy J dy 1 a ab (x, x') , (52) 

C a C b 

where C a and Cb are the cross-sections of the leads and 
a a b is the non-local conductivity for graphene^Ir— , 

Vab(x, x') = ^Tr [a a G(x, x')* b G*{x, x')] . (53) 

Here, the vectors a and b are unit vectors in the di- 
rection of the corresponding lead and point into the in- 
terior of the system (cf. Fig. [5]). In App.O we show 
that this method to compute g a b is equivalent to an ap- 
proach based on a Fisher-Lee type formula 4 ^ generalized 
to graphene, which we also derive in App.lCl 

Inserting the semiclassical Green's function © into 
Eq. (f53")) . we obtain a double sum over classical orbits 
7 and 7' that lead from point x' in lead b to point x in 
lead a: 

a ab (x,x') « (^j'EW/^"^ 
77' 

(54) 

with SSj t Y = hk E {L 1 — Ly), <5/i 7l7 ' = /z 7 — /Lty and 

F 7i7 , = Tr (a a K y a b K\?) . (55) 

We note that the entire effect of the graphene pseudospin 
and valley dynamics is included in the factor Y ltl '. The 
main difference of this factor from Z ltJ i in Eq. ([T7]) is 
that the pseudospin propagators have to be multiplied 
before tracing. 



A. The average quantum conductance 

1. Quantum chaotic Schrodinger billiards 

First we focus on the average conductance, 
{gab)k E = (gab) where as in Sec. lIIII the averaging 
is performed over an energy window that is classi- 
cally small but still contains many quantum levels. 
Based on semiclassical approaches for the Schrodinger 
cas o 24 ' 27 i 41 ~ — and random matrix theor y 44 ' 45 coinciding 
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universal predictions for the average quantum conduc- 
tance were made for chaotic cavities for the unitary, 
orthogonal, and symplectic universality class: 



(dab] 

9o 



M n M h 



M a + M b - 1 + 2//3 



(56) 



Here (3 depends on the universality class and the consid- 
ered random matrix ensemble, respectively, namely f3 = 1 
for the circular orthogonal ensemble (COE), (3 = 2 for the 
circular unitary ensemble (CUE) and (3 = 4 for the circu- 
lar symplectic ensemble (CSE). The conductance quan- 
tum is defined as go = e 2 /h and M a / b is the number 
of propagating channels in the corresponding lead. For 
large M a / b , we can expand (g a b) as 



{g a b) _ M a M b 



M a + M b 



M a M b Q ( M -r 
(3 J {M a + M b Y +U V«/ b ) 

(57) 

Semiclassically, just as for the spectral form factor, rele- 
vant contributions to the average conductance (|57|) arc 
due to orbit pairs that have a small or vanishing ac- 
tion difference 8S lri '. In fact for generic pairs with 
<55 7 . 7 ' ~ hk E L the fastly oscillating exponential in (|54[) 
leads to zero average conductance. Here, similar to the 
case of F(t), we compute contributions to (g ab ) due to 
diagonal and leading off-diagonal terms in Eq. (|54|) in the 
next two subsections. 



2. Diagonal contribution 



We begin with the diagonal contribution, where 7' = 7 



and 5S. 



7.7 



7.7 



0. We note that the bound- 
ary matrices in Eq. ([5]) do not mix valley and sublatticc 
pscudospin. Thus we can separate the graphene pseu- 
dospin propagator into a sublattice part and a unitary 
valley part K^, which we trace out immediately giving 
Tr (K^K^} = 2. The pseudospin contribution in diago- 
nal approximation then becomes 



Y 7 = Tr (a a K s 7 a b K^) Tr (iqK?) 
= 2Tr (a a K s j( j b K^) . 



(58) 



The sublattice part depends on the angles of the incoming 
and outgoing trajectory segments (cf. Figj5]) leading to 



icos(6> a )cos(0 b ) , 



(59) 



as shown in App. lA3l This result for the pscudospin 
trace within diagonal approximation holds irrespective 
of the specific edge types involved in the various reflec- 
tions along 7 (in contrast to the diagonal contributions 
to the spectral form factor). We now compute the full 
diagonal contribution to Eq. (|54[) as single sum over or- 
bits with fixed outgoing and incoming angles employing 
the classical sum rule^i 



£ \d,\ 2 s(t-t,) = 



v F T.(k E ) 



d6 a d9 b e- T/Td (60) 



to transform the sum into a triple integral over the time 
the particle spends in the cavity and entrance and exit 
angles. Here T<j is the dwell time, i. e. the time a classical 
particle typically spends within the cavity. In terms of 
the cavity area and the lead widths W a and W b , the 
corresponding escape rate is approximately given by (see 
Ref.|36| and references therein) 



1 



ttA 



(W a + W b ) 



(61) 



Additionally we find for the energy surface £ in Eq. ((60)) 

£(M = J d 2 x J d 2 kS(k E -k) = 2tt Ak E . (62) 

Combining Eqs. ([541 159l6"2")) we get for the diagonal non- 
local conductivity 



(a D (x,x')) 



" \4wh) 

hk E 

v p 2-kA 
2k P 



t/2 tt/2 

dd a J d6 b 8cos(9 a )cos(6 b ) 

-TT/2 -7T/2 



dTe- T ' T - 



9o 



7T W a + W b 



(63) 



(64) 



After integration over the lead openings in Eq. (|52"|) we 
obtain 



{9D) 
.90 



2k E W a W b 



M a M b 



tt W a + W b M a + M b 



(65) 



without edge signatures and in agreement with the lead- 
ing order term in expression (|57l) . In Eq. (|65[) we have 
approximated the numbers of propagating lead channels 
M a and M b by 



M, 



2k r 



a/b 



■w. 



a/b 



(66) 



accounting for the two valleys. For identical leads (M a = 
M b = M) we have (go) /go = M/2, half of the maximum 
possible conductance, as expected for a classical particle 
being randomly scattered. 



3. Loop contributions: Weak localization 

Wc now evaluate the leading order quantum correction 
to the average conductance, namely the loop contribu- 
tions that are responsible for the weak localization (WL) 
correction in (g ab ) ■ The corresponding orbit pairs (7, 7x ), 
sketched in Fig.[5j are similar to the ones considered in 
Sec. IIII1 however here the orbits are not periodic but they 
are open orbits connecting one lead to the other. Wc di- 
vide the orbits into the legs l\ and I2 where 7 and 7 X 
both have the same propagation direction and the loop 
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FIG. 5. Sketch of a loop pair that contributes to the weak 
localization correction to the average conductance^. One of 
the orbits crosses itself with an angle e, the other orbit avoids 
this crossing but apart from this follows the first one expo- 
nentially close. 



part L, where the propagation directions are opposite 
(cf. Fig. [5]), and label all quantities correspondingly. For 
vanishing Ehrcnfcst time, the effects of a finite encounter 
time can be neglected here. In App. lA4l we compute the 
trace of the pseudospin matrices of the loop pairs and 
find 



Y, 



7,7x 



!C0s(# a ) COs(#b) 



(67) 



-1 for odd N aCyL , 

(-1) Nl cos{AKk L + 2tf L ) for even N ac , L . 

For the summation in Eg. (|54l) . Y 7nx has to be aver- 
aged according to the edge characteristics of the sys- 
tem. We make use of the crgodic dynamics of classically 
chaotic cavities, particularly the fact that long orbits hit 
all points on the boundary with the same probability. 
Therefore we average Ky i7x over the boundary conditions 
corresponding to the different edge types. Here again the 
ac edges play a crucial role: Only classical orbits that hit 
ac edges an odd number of times during the loop part 
lead to a finite contribution to the average conductance. 
For pairs with odd A^ aCj i we have an average pseudospin 



contribution (V 7 



-8 cos( 



as follows di- 



rectly from Eq. (|67|) . On the other hand for even N ac ^ 
we obtain (Yy J y x ) = O 4 ^. To see this we first assume that 
all reflections occur at ac edges, i.e. Nl = N aCi L and 
&l = 0. Since for a generic boundary KAl = mr/3 [cf. 
Eq. (|54[) ] with n a random integer, the cosine in Eq. (|67[) 
is zero on average. The result of this averaging is not al- 
tered by the inclusion of reflections from zz and im edges. 
Using the sum rule (|60|) and replacing the 7 sum by inte- 
grations over T, 9 ai 6 bl and the crossing angles e, we find 
then for the loop correction to the non-local conductivity 



(a L (x,x')) = 



2c V]?kft , 
it 3 HA 



7T OO 



xP(s,T)e 



Re / de / dTe-V- T —^' T - 



(68) 



2T min (e) 



Note that taking the real part and including an addi- 
tional factor of two accounts for the fact that every orbit 



pair in the sum has a partner pair, where 7 and 7 X are 
interchanged leading to the complex conjugated term. In 
Eq. (f6"8j) P denotes the density of loop pairs that hit ac 
edges an odd number of times along the loop L. A weak 
magnetic field that causes a flux $ through the area of 
the cavity leads to dephasing on a timescale 



$2 



(69) 



with the magnetic flux quantum $0 = h/e and a system 
specific parameter sec Eq. (|B3|) . In App.lBl we show 
how this affects P. The modified loop density is then 
given by 



P(s,T) 



9 

2ttA 



T-T min ( £ ) 



dT L [T~T L - T rain (e)] (70) 



T min ( e ) 



x e 



-[Tz,-T min (e)]/T B 



1 - e 



-[T L -T min (e)]/T a , 



where T ac is the typical time between two reflections from 
intervalley scattering ac edges [see Eq. ([28)1 ] and we used 
the probability for an odd number of ac reflections (|49p . 
Further T m ; n is the minimum time to form a closed loop 
[see Eq. dM])]. Evaluating the integrals in Eq. (|6"8"]1 in the 
limit AT d > ln(w F fc E /A), we find2i^ 



{9L) 
90 



M a M b 



1 



1 



{M a + M b f 



\ _|_ Ik 



1 



(71) 



This central result represents the leading order quantum 
correction in inverse channel number (known as weak 
localization correction) to the average magnetoconduc- 
tance of a ballistic grapheme cavity. Three competing 
timcscalcs govern (gi) : the armchair scattering time T ac , 
Eq. (|28|) . the dwell time Td, Eq. (|B"Tj) . and the magnetic 
time Tb, Eq. (|69| . Equation (|7T|) describes the WL - no 
WL crossover as a function of the amount of armchair 
scattering and the magnetic field. Unlike for the spec- 
tral form factor, Eq. (|5T|) . there is no dependence on the 
Fermi momentum k E (for fixed numbers of propagating 
lead channels). For vanishing $ the size of the WL is 
given by 



r (9l) 
lim 



M„.M h 



1 



(72) 



(M a + M b ) 2 l + T ac /(2T d ) ' 
describing the T ac dependence of cf. inset in Fig. [5] 

a. ac billiards - complete valley mixing 

For very strong armchair scattering, T ac <C Td, the 
well known Lorentzian magnetoconductance profil e) 24 ' 41 
is restored, 



lim M 



M„M h 



1 



(M„ + M b y 1 + T d /T B ' 



(73) 



describing the usual GOE — > GUE transition, where the 
total WL signal has a size of approximately —1/4 for 
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identical leads in the limit of large channel numbers. 
Since the valleys are completely mixed, the TRS T y is 
active at $ = 0, leading to the orthogonal universality 
class. The magnetic field breaks T y , driving the system 
into the unitary class. 



magnetic dephasing time, the semiclassical calculation 5 ^ 
is unchanged. Thus performing the integrals in Eq. (|68[) 
for finite Ts/Td, the WL is suppressed in the same way 
also here, namely in Eq. (|7ip we get an overall exponen- 
tial suppression 



b. Weak ac scattering 

For very weak ac scattering, T ac ^> Eq. (f7Tj) be- 
comes 



{9i 



lim 



M„Mh 



1 



2T H 



(M a + M h y (1 + T d /T B ) 2 T ai 



(74) 



that is the WL correction vanishes linearly with 2d/T ac , 
reflecting (partial) effective TRS breaking. The overall 
TRS associated with the antiunitary operator T y is in- 
tact for zero magnetic field, implying orthogonal sym- 
metry. However, in the limit of vanishing ac scattering, 
Id/Tac — > 0, the Hamiltonian is block diagonal and T y 
is ineffective since it connects states from different val- 
leys. Within each valley the effective TRS associated 
with 7o is broken by the zz and im boundary condi- 
tions thus (<7l) is suppressed but never becomes positive. 
This is in contrast to the case of disordered bulk samples 
with small intervalley scattering rate, where the symme- 
try corresponding to 7o is not broken leading to weak 
anti-localization^. That means we have two (nearly) de- 
coupled unitary subsystems that lead to vanishing (or 
strongly suppressed) WL. 

We point out that the magnetoconductance curves 
described by Eq. (|74[) are squared Lorentzians and thus 
steeper and more narrow than Lorentzians with the 
same height (cf. dashed blue line in Fig-EJ) • That means 
the magnetic dephasing is essentially assisted by the 
partial TRS breaking. 

In Fig. [5] we summarize our findings graphically. The 
main panel shows the dependence of the WL correction 
on the magnetic flux, or y/T d /T B = y/T<i/C$/®o, 
respectively, for different values of the ratio 7d/T ac . 
As the ac scattering becomes stronger, the magne- 
toconductance gets closer to the Lorentzian behavior 
known from usual 2DEGs. The inset shows the crossover 
as a function of T&/T ac for the case of zero magnetic flux. 

To obtain Eq. ([TTj) . we have assumed that the Ehrcnfcst 
time 



T E = i h\(v F k E /X) 



(75) 



is much smaller than the dwell time Td and thus neglected 
terms of order Ts/Td or higher. Te is the time scale af- 
ter which a minimal wave packet in a chaotic cavity has 
spread over the whole system, so that quantum interfer- 
ence is possible^. For Schrodinger systems it is known52, 
that a finite Ehrenfest time leads to an exponential sup- 
pression of the WL correction. Since for graphene the ac 
scattering enters the integral (|70l) in a similar way as the 



K9L) 
9o 



M a M b 



(m„ + M b y 



-T E /T d 

TUT 



-T E /T d 



Ik 



(76) 



4- Comparison with numerics 

In Fig. [7] we compare our analytical results with the 
conductance obtained from numerical tight binding sim- 
ulations. Here we consider a graphene cavity with the 
shape of a desymmetrized half stadium giving rise to 
chaotic scatterin g 52 ! 53 . We model infinite mass bound- 
aries by adding a smooth mass term close to the edges 
of the system (for details see Ref. fl9l ). In order to tunc 
T ac we cut a hole into the smooth mass so that locally 
a small region with armchair edges is present. The sys- 
tem is shown schematically as an inset in Fig.[7J We 
then calculate the magnetoconductance numerically (the 
inset shows three examples) using an adaptive recursive 
Green's function method 5 ^ and obtain (gr) as the differ- 
ence of (#)($) at zero flux and $wl.4 - 1.7$o- The 
main panel in Fig. [7] compares the absolute value of the 
WL correction obtained in this way with the semiclassical 
prediction (|72|) . We find that the numerical data repro- 
duces our semiclassical theory very well (dashed line). 
The agreement is even better if we assume an effective 
offset in W ac of five lattice constants, AW ac = 5a (solid 
line). This offset can be explained by the fact that our 
smooth mass edges cannot completely avoid intervalley 
scattering from armchair edges or sharp edges at the 
boundary, e. g. at the lead mouths. We assign further 
deviations of the numerical from the anayltical results to 
the limitations of the Dirac equation and the semiclassi- 
cal approximation. 



B. Universal Conductance Fluctuations 

So far we have considered quantum effects in the aver- 
aged conductance, i. e. the smooth part of g. Now we ad- 
dress quantum fluctuations around this average value, of- 
ten denoted as universal conductance fluctuations (UCF) 
because their variance does not scale with the average 
conductance of a disordered or chaotic system. The UCF 
in diffusive graphene samples have been investigated in 
several theoretica l 55 ' 56 and exp er iment a l 1 ' 5 7 ~ — studies. 
Here we consider the case of ballistic graphene cavities. 
To quantify the size of these fluctuations, we calculate 
the variance of the conductance g ab which is identical to 
the covariance of the resistances g aa and gbb- 



var(g ah ) = (g aa 9bb) - (9aa)(g 



(77) 
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FIG. 6. Absolute value of the WL correction, — (<?l), as 
a function of the magnetic flux $ through the system, ob- 
tained from our semiclassical theory for identical leads. The 
full lines correspond to different values of the effective arm- 
chair scattering strength Td/T ac (from bottom to top) : 0.1 
(black), 0.3 (red), 0.5 (blue), 1.5 (turquoise), 5 (orange) and 
oo (green). The dashed blue line shows a Lorentzian curve 
~ [l + ($ 2 T d )/($oC)] _1 f° r comparison. The inset shows the 
crossover as a function of Td/T ac at zero magnetic flux ac- 
cording to Eq. (f?2"j) . 
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FIG. 7. Comparison of the semiclassical theory for weak lo- 
calization with numerical tight binding simulations. For the 
numerics we use a system (with leads W a = Wt = 40 a) 
with a smooth mass edge modeling infinite mass boundaries^ 
and open this mass boundary in order to tune W ac - Fitting 
Eq. (|71[) to the numerically calculated magnetoconductance, 
we obtain the total height of the weak localization signal 
[1 + Tac/(2T d )] -1 as a function of W^. The inset shows 
(from bottom to top) three examples for W ac = (green) , 
50a(black) and HOa(red) . We find £ = 0.17T d and for the 
limiting value (<7i)($ = 0, T ac = 0) = 0.27 go from a fit to the 
numerical results for a system without a mass boundary that 
shows very strong intervalley scattering. The solid (dashed) 
line in the main panel displays the crossover as a function of 
Wac obtained from the semiclassical theory, Eq. Q72p. includ- 
ing an offset Aiy ac = 5a (AVKac = 0). The error bars are 
centered around the numerical data. 




FIG. 8. Schematic view of orbit quadruplets that give rise 
to the conductance fluctuations in leading order, see Ref.l60l. 
Here we have two (avoided) crossings of orbits from different 
leads. The pieces of the orbits before, after and between the 
crossings are labeled with letters ntor. 



Here the semiclassica l 53 ! 60 and RM T 44 ' 45 ' 61 prediction for 
many propagating lead channels reads 

var(g ab ) 2M a M b 



9l 



/?(M„ + M b ) 4 



(78) 



for a scalar Schrodinger equation. We consider the first 
term in Eq. (f77|) . the averaged product of resistances. Ac- 
cording to Eq. (|54[) . multiplication of two resistivities re- 
sults in the fourfold sum over classical orbits 



2 2 



77' (a) pp'(b) 



16tT 2 H 2 j 

x exp [idS^ tl >/H + iSS PiP ' /h + iSp 7n >ir/2 + i5p p . P 'iT /2] . 
with 



X. 



i,i',p,p' 



Tr 



iKj<j 



Kl, ) Tr 



(79) 



Note that the orbits 7 and 7' originate from the same 
point x' a in lead a and end at the same point x a also 
in lead a. The orbits p and p' on the other hand both 
begin at x' b in lead b and end at x b in lead b. Here the 
averaging procedure gives only a significant result when 
the sum of the action differences is small or vanishing. 
One possibility is that the action differences of the orbits 
from one lead, i. e. SS ln ' and 5S PyP ', are both small indi- 
vidually and the orbits from different leads are uncorre- 
cted. These contributions are fully contained in the sec- 
ond term, namely the product of averages, and thus they 
cancel out. Therefore we have to consider orbit quadru- 
plets in which only the complete phase is small, but not 
the individual phases. Neglecting effects due to a finite 
Ehrcnfcst time, the dominant pairings arc those sketched 
in Fig. [51 as shown by Brouwer and Rahav^S. Two orbits 
from different leads, say 7 and p, cross each other with 
a crossing angle £\ and then cross each other again, now 
with a crossing angle £2- The two other orbits follow the 
first ones exponentially closely, but avoid the crossings. 
For instance 7' follows 7 to the first crossing, then it fol- 
lows p to the second crossing and finally again 7 back to 
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the lead [Fig.[8]a) and b)]. We label the individual pieces 
with to, n, o,p, q, r as marked in Fig.[5]a) and denote the 
combination of the pieces r and q~ 1 , which is q traversed 
in opposite direction, by the loop part L. Note that also 
quadruplets where two orbits from different leads cross 
first and then avoid the second crossing (or vice versa) 
while the two other orbits avoid the first crossing and 
cross then (or vice versa) give rise to a contribution to 
the UCF in leading order [Fig.|8]c) and d)]. 

For the details of the calculation of X~ (n ^ PtP ' we refer 
to App. lA 5l We note however that, in contrast to the 
individual traces in Eg ([75)) . their product X in ^ p ^ p i de- 
pends only on the reflections along the loop part L. In 
fact we find X 



7,7 ,P,P 



if the total number of reflec- 
tions from ac edges during the loop part is odd, and for 
an even number of ac reflections we get 

X-y,y,p,p' = 64cos(6> a ) cos(8' a ) cos(8 b ) cos(8 b ) 

xcos 2 (2KA l + $ l + N l tt/2). (80) 

In order to perform the orbit sums, we make use of the 
fact that the two self-intersections are independent and 
thus their contributions factoriz e 42 ' 60 so that we get 



var(cr Q 6 



2e 2 v F k E W a W b 



TT 3 HA 

oo 



x Re j dei 

2T„ 



dTi e 



(81) 



-[T!-T mill ( £l )]/T dp i5S( El )/fi 



n(ei) 
oo 



x Re J de 2 J dT 2 e~ [T2 

2T min ( £2 ) 

x P ac (e 1 ,E 2 ,T 1 ,T 2 ). 



,(e 2 )]/T ie i5S(e 2 )/h 



Here P ac (£i) £2, Ti, T 2 ) is the density of two uncorrected 
crossings such that armchair pieces are hit either an even 
number of times during both the q and the r part of the 
loop or an odd number of times during both parts 

Ti-T min (£i) T 3 -T min (e 2 ) 



P„c(ei,e2,Ti,r a )M -z- 
1 nA 



dT„ 



dT r 



x [7\ -T q - T„ 



T mi „( £l ) 

i(ei)][T 2 -T r -T a 



T min (£ 2 ) 



x P a e c [T f; + T r - T mi „(£i) - T min (£ 2 )] 



(82) 



To obtain Eq. (|8"Tj) we have used the sum rule (jSTTf twice, 
included the averaged pseudospin traces {X^.^i . p . p >) for 
the contributing quadruplets and integrated over the in- 
coming and outgoing angles as well as over the lead open- 
ings. Furthermore we included a factor of two due to 
time reversed quadruplets, i. e. such quadruplets, where 
the orbits from different leads (e. g. 7 and p') traverse the 
corresponding part of the loop (q) in opposite direction 
rather than in parallel [panel b) vs. a) and d) vs. c) in 
Fig.|8]. For zero magnetic field, both contribute equally, 
and analog as for Eq. (fTTj) we get for the variance of the 



conductanc o 24 i 48 



v&r(g ab 



2M 2 M 2 



9l 



1 



(1 + 2T d /r ac r 



,(83) 



which is the central result of this section. The crossover 
described by this equation is different from the one for 
the average conductance, Eq. ([72]) . Also here we have no 
explicit dependence on k E , but as for the spectral form 
factor, degeneracies have to be taken into account when 
interpreting the result. For identical leads we have a 
crossover from 1/4 to 1/8 for an increasing probability of 
armchair scattering. This agrees with a transition from 
two degenerate unitary subsystems (4xCUE) to one or- 
thogonal system (COE). 

For a finite magnetic flux through the cavity we have 
to modify the densities accordingly as in the previous 
sections. However, not all quadruplets are affected by 
the magnetic field, because there is an additional phase 
due to the flux only if the loop part is traversed in oppo- 
site directions by the different orbits. Thus we get one 
half of the zero flux result plus the contribution from the 
quadruplets that are affected by the magnetic dcphasing 
for the variance: 



var(5 Qh ) 



M 2 M 2 



9l 



~ (M a + M b y . y (l + i Td /T B + 2j T d /T ac ) 2 

(84) 

Even though all orbit quadruplets arc affected by the 'de- 
phasing' due to the armchair scattering, the magnetic 
flux has an effect only on one half of the orbit quadru- 
plets and both mechanisms lead to a similar behavior of 
the variance in the final result (|84[) . For identical leads 
and T ac —> 00, Eq. (|8~4"|) describes a crossover from 1/4 
to 1/8 with increasing flux $, i.e. from two degenerate 
unitary subsystems (4xCUE) to two independent uni- 
tary subsystems (2 x CUE) , and for T ac — > from 1 /8 to 
1/16, i. e. from an orthogonal (COE) to a unitary (CUE) 
system, in agreement with our findings for the average 
conductance and the spectral form factor. In Fig. [9] we 
depict our semiclassical result (f8~4")l . We show the mag- 
nctodcpcndence of the conductance variance for differ- 
ent values of the ac scattering time T ac (main panel) as 
well as the dependence of var(g) on the ratio T^/T ac for 
zero (inset, upper black curve) and finite (inset, lower 
blue curve) magnetic flux. The crossover as a function 
of both the ac scattering and the magnetic flux is here 
faster than for the WL correction to the average conduc- 
tance, because the denominators in Eq. (|84p are squared 
as opposed to those in Eq. (fTTT) . 

Fig-Hni compares the semiclassical result with numer- 
ical tight binding calculations. The system used for the 
numerical study is the same as we used for the average 
conductance. We find that Eq. (|84p provides a good ap- 
proximation to the numerical data. 
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FIG. 9. Variance of the conductance as a function of the 
magnetic flux $ through the system obtained from our semi- 
classical theory for identical leads. The curves correspond 
to different values of the relative armchair scattering strength 
Td/T a c (from top to bottom): 0.0 (black), 0.1 (red), 0.3 (blue), 
0.5 (turquoise), 1.5 (orange) and 5 (green). The inset shows 
the crossover as a function of Td/T ac at zero magnetic flux 
$ = (upper black curve) and at finite flux $ = 5 ^o/a/^VC 
(lower blue curve) [cf. Eq. (J83J) and Eq. d84}] . 
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FIG. 10. Comparison of our semiclassical theory, Eq. 
with numerical tight binding simulations. For the numerics 
we used the same system as for the average conductance. The 
main panel shows the crossover as a function of VK a c obtained 
from the semiclassical theory Eq. H84|) (solid line) and from 
the numerical data (lines with error bars) at zero magnetic 
flux $ = (upper black curve) and at finite flux $ = 1.4 "I>o 
(lower blue curve). The inset shows (from top to bottom) 
the dependence of the numerical data on the magnetic flux 
(lines with error bars) for systems with W ac — (green), 
50a (black) and 110a (red). We find that £ « 0.25 T d fits 
the data well. For the solid lines at W ac = and 50 a, we 
added a constant offset of 0.011 to the result of Eq. Q84p , the 
dotted lines show the result without any offset. 



C. Shot noise 

As another important transport property, we briefly 
discuss the shot noise of a ballistic graphene cavity. For 
graphene ribbons with rather large width to length ra- 
tio, where transport at very low energies may be domi- 
nated by evanescent modes, shot noise has been studied 
theoreticall y 62 ' 63 and experimentall y 64 ' 65 . We focus here 
on the case of chaotic scattering at larger energies, ap- 
plying our semiclassical theory. Usually the shot noise is 
quantified via the Fano factor J-, defined as the ratio of 
the shot noise power S of the system and the Poissonian 
shot noise power Sp corresponding to charge carriers be- 
ing transmitted in an uncorrelated way from one lead 
to the other. In terms of the transmission amplitudes 
t = tab the Fano factor, that describes the quantum sup- 
pression of the shot noise with respect to the classical 
value, is given bj^ 6 - 



T 



S_ 



Tr (ii* - ttHtf) 



Tr(ttt) 

Using RMT for the scattering matrix one obtains^ 5 - 

M a M b 



T 



(85) 



(86) 



In App.0 where we discuss the generalized Fisher-Lee 
relations, we show that 

Tr (tit) = ~ J dy J dy 1 a ab (x,x') , (87) 

c a c b 

Tr (tthtf) = Jdyi jdy 2 Jdy[ Jdy' 2 a ab (x 1 ,x 2 , x\,x' 2 



C a C a Cb Cb 



with 



a ab (x 1 ,x 2 ,x' 1 ,x' 2 ) (89) 
= Tr[a a G(x 1 , x' 1 )a b G\x 2 , x' 1 )a a G(x 2 , x' 2 )a b G\xi, x 2 )] . 

Inserting the semiclassical Green's function <j2j) we get 
the fourfold sum 

a ab (x 1 ,x 2 ,x' 1 ,x' 2 ) = (j^j (90) 

x Y,Y, Tl (eaK^o-bK^c-aKpO-bKl^ D 7 D Y D P D P , 

77' PP' 

x exp [iSS 7i7 '/H + iSS PtP '/h + i<5/i 7!7 '7r/2 + iS/ip^'ir /2] , 

similar as for the UCF. However, for the UCF the con- 
tributing classical orbits came in pairs with the same 
starting and end points. Here, the orbits with the same 
starting point have in general different end points and 
vice versa. In words, the trajectory 7 starts at a point 
x\ in lead b and ends at a point xi in lead a, the tra- 
jectory 7' starts also at x[ but ends at a point x 2 in 
lead a. The trajectories p and p' both begin at a point 
x' 2 in lead b and p ends at x 2 while p' ends at x\. The 
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FIG. 11. Sketch of the orbit quadruplets contributing to the 
shot noise in leading order. As opposed to the UCF only one 
crossing is present. Orbits that begin at the same point, end 
in general at different points. 



diagonal contribution corresponds to quadruplets where 
either 7 = 7' and p = p 1 or 7 = p' and p = 7'. This is 
only possible if either the end points are equal x\ = x-2 
or the starting points x[ = x' 2 . One can show that in 
this case Tr (ftt) = Tr (ttUtf)£L§2., and the Fano factor 
T is zero in diagonal approximation. This is the classical 
limit, where no quantum shot noise exists. The leading 
order contribution to T is known to stem from diagrams 
as sketched in Fig.QjJZr— . For these quadruplets the val- 
ley part of the pseudospin trace does not depend on the 
edge type at all, since 

Tr (k;k;1k;K v J) = Tr {K v q K v K v jK^K v r K v p K^K^) 



= 2. 



(91) 



Also for the sublattice part, we do not find such an edge 
dependence. A longer calculation similar to the one for 
the average conductance and the UCF results in 



Tr ^ a K s , ( a b K;\a a K s p a b K p , 
= 32 cos(0i ) 003(61;) cos(6> 2 ) cos(6» 2 ) . 



(92) 



Hence to leading order in the inverse channel number 
the Fano factor of chaotic graphene flakes does not de- 
pend on the edge characteristics. The orbit sums can 
be performed exactly as in the case of usual Schrodingcr 
systems giving 



Tr (it 1 - tthtf) 



(93) 



which, together with Eq. (|55|) . yields the leading order 
term of Eq. (|5B")) . To conclude, at least to leading order 
in inverse channel number, the shot noise does not exhibit 
any imprints of the edges of a chaotic graphene cavity. 



V. CONCLUSIONS 

The major part of theoretical work on quantum in- 
terference phenomena in mesoscopic graphene systems 
has focused on the diffusive regime. However due to the 
permanent progress in manufacturing techniques it be- 
comes possible to fabricate clean structures, where the 



main scattering mechanism is due to reflection from the 
system's edges. In this work, we provide a comprehen- 
sive semiclassical study of quantum interference effects 
in ballistic graphene nanostructures, particularly focus- 
ing on the effects of the edge characteristics. To this end 
we have started from a general analytic expression for 
the full Green's function of a graphene cavity with arbi- 
trary edge structure^ and generalized semiclassical meth- 
ods that have been originally developed for Schrodingcr 
systems to graphene nanostructures. 

As a first application we have have studied the fluc- 
tuations of the density of states of closed graphene bil- 
liards with chaotic classical dynamics. Since the DOS 
of such billiards is very complex, we have considered its 
statistical properties, namely the two-point correlations 
of the DOS via the spectral form factor F(t). Already in 
the diagonal approximation we find a strong dependence 
of F on the specific edge structure. We find that the 
total amount of intervalley scattering armchair bound- 
aries sets the timescale for a crossover in the form factor, 
where the asymptotic limits agree with the RMT predic- 
tions for two degenerate unitary subsystems (4xGUE) in 
the absence of intervalley coupling and a single orthog- 
onal system (GOE) in the limit of complete valley mix- 
ing respectively. We also have derived the explicit edge 
dependence of the simplest off-diagonal corrections, con- 
firming our interpretations of the diagonal contributions. 
Even for a significant amount of intervalley scatting ac 
edges, we predict that the small t correlations deviate 
significantly from the GOE result. Our main results for 
the DOS, Eqs. (|30| and (|5Tj) . thus suggest that in exper- 
iments, the effects of effectively broken TRS should be 
visible in the correlations of the DOS on a scale of many 
mean level spacings rather than in the usually consid- 
ered next nearest neighbor statistics, where intervalley 
coupling is essentially too effective to see those effects^. 
Our results do not support the experimental indication of 
GUE statistics in the nearest neighbor level distribution 
in Ref.0. The latter would suggest extremely weak inter- 
valley scattering or some other, yet unknown, symmetry 
breaking mechanism. 

In the second part of the paper, we have further derived 
analytical expressions for the most important transport 
properties of open ballistic graphene cavities in the semi- 
classical limit, such as the average conductance, the uni- 
versal conductance fluctuations and the Fano factor. We 
have calculated the classical value (diagonal approxima- 
tion) and the leading order quantum corrections, i. c. the 
weak localization correction, for the average conduc- 
tance. While the classical conductance does not depend 
on the edge structure of the cavity, the WL correction is 
again governed by the effective armchair scattering time. 
Equation (fTTj) gives the WL correction of a graphene flake 
in terms of its microscopic edge properties. The WL is 
suppressed if intervalley coupling is absent, i. e. there is 
no weak antilocalization in this case, because the effec- 
tive intravalley TRS is not preserved in the presence of 
edges. With increasing ac scattering strength, the inter- 
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valley TRS becomes effective again and the WL correc- 
tion reaches the COE value for very strong intervalley 
coupling. For the UCF we find corresponding crossovers 
between the unitary and the orthogonal random matrix 
ensembles in agreement with our findings for the WL and 
the spectral form factor. Finally we have shown that the 
shot noise power and thus the Fano factor are not affected 
by the edge structure to leading order. Our main results 
for the transport properties, Eqs. ([71]) and (|84|) explain 
part of our earlier numerical worki^, and provide a theo- 
retical footing for experimental and numerical studies of 
graphene samples in the ballistic regime. 
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Appendix A: Calculation of pseudospin traces 

1. Pseudospin trace for the form factor F(t) 

We calculate the trace of K 1 , Eq. (fT3|) , for a given clas- 
sical orbit 7 taking into account reflections from edges of 
infinite mass, armchair, and zigzag type. We separate 
the pseudospin propagator into a sublattice and a valley 
part, K y = ® K^, and consider first the valley part 
K^. It is traceless if the number of ac reflections is odd, 
and otherwise we have [cf. Eq. ([5])] 

k; = ±7 jv(2*a,+*,)t, (ai) 

and thus 

TrK^ = ±i N 2 cos(2A'A 7 + 7 - Att/2) , (A2) 

with A 7 and $ 7 as defined in Eqs. (|T5|) and (jTij)) . The 
sign depends on the exact sequence of zz, ac, and im 
reflections and the individual signs in Eq. ([5]) . Every pair 
of successive ac reflections gives rise to 

T v e 2lKx *T y e 2lKx ^ = e 2^(z.-i-^) (A3) 
and thus to A 7 . 

For the sublattice part, we begin assuming that only one 
reflection along 7 is at a zz edge, without loss of gener- 
ality the very first one, and afterwards an even number 
of reflections occur at ac edges and the rest at im edges. 
Then we can cast in the form 

Kl = i N a^ 1 e^- e ^a tl (l + a ai , a „) (A4) 



and therefore 

TrA' 7 = 2 l 2A '- 1 cos [9 7 —9i+(p ai , aN —<pt 1 — {N—l)ir/2] , 

(A5) 

where tp v is the polar angle of the vector v. Geometric 
considerations give further 

fa 1 ,a N — <Pt 1 =0i—it/2 mod 27T (A6) 

and thus we get 

TrA 7 = 2i 2N - x cos (6> 7 - Att/2) . (A7) 

Similarly we write for the case of two zz reflections 

Kl l = l N af\^ ) °*a t2 af ) e< 1) °--a tl {l + a ai . aN ) 

(A8) 

leading to 

TrA 7 = 2i 2N+2 (-l) Nil)+1 

x cos (0W - 0(U + tp ti - Vt2 - Att/2) .(A9) 

Once more employing a geometric relation 

e^-e^+(p tl -cp t2 =0 7 -(N (1) + l)7r mod27r (A10) 

we find 

TrA 7 = 2i 2N ~ 2 cos (6> 7 - Att/2) . (All) 

The treatment of these two special cases can be combined 
to find the pseudospin trace for an arbitrary number of 
zz reflections. Absorbing all phases in the prefactors into 
/ 7 , we finally obtain Eq. (fT3")) . 

2. Trace Z 7 , 7x for the loop contribution to F(t) 

From the discussion above, we know already that 
TrA 7 = TrA' 7x = if the total number of armchair 
reflections is odd, so we can concentrate on the opposite 
case. By using Eq. (|3"2"|) . we find 

cos(0 7 + Att/2) cos(0 7x + Nn/2) 
= cos 2 (0 L + Att/2) cos 2 {0 r ) - sin 2 (0 L + Att/2) sin 2 (0 fl ) 
= (-1) Nr+1 + 0(s 2 ), (A12) 

where we made use of the expressions equivalent to 
Eq.fZB 

0l = (N+, L - A_ :L + 1)| - w l tt + O(s) , (A13) 
Or = (N+, R - A_ ifl - 1)| - Wr tt + O(e) . (A14) 

Since in the semiclassical limit the main contribution to 
Fl comes from pairs with small angles e— , we keep only 
the leading term (-1)^+1. With that Eq. ([33]) for the 
loop pairs with even A^ ac follows from Eq. (|13[) . 
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3. Trace V 7 for the diagonal contributions to the 
conductance 



For the sublattice part of Y~ ( we consider orbits that 
hit armchair and infinite mass edges. Each reflection ef- 
fectively contributes a term a z e iC7 * and thus 



Tr (a a K y a b K^) 
= 21t[a afe <fl ^.(i + w ) ol (i + w ) e -^. ^] 
= (-l) w 2Tr [e- 2l6 ^{o a o b + a a o aiX ,o b a aiX ,)\ 
= (-l) w 8 cos(26» 7 + if a - tp aiX ')cos{(p b - ip aiX ') ,(A15) 

where ip a , ip b , and f aiX ' arc the polar angles of the vec- 
tors a, b, and (ati — x') respectively. We define the angles 
of the incoming (6b) and the outgoing (6 a ) segment of the 
trajectory 7 relative to the lead orientation (cf. Fig. [5]), 



1a = <£xa N — (fa + K ■ 



(A16) 

(A17) 



and make use of the relation 

26 1 = Ntt + ip aix , - <p xaN mod 2tt . (A18) 
so that we finally obtain 

Tr (a a K 7 a b K\) = -8cos(0 a ) cos(6 b ) . (A19) 

By using the same geometrical relations as in App. lA ll 
one shows that this holds also if reflections from zigzag 
edges are involved. 



Trace Y~. 



for the weak localization correction 



For the sum of reflection angles as defined in Eq. (fT4")) 
wc have 



#7 — Oiy + 6l + 6[ 2 , 
^7x = @h ~ ®l + 6i 2 . 



(A20) 
(A21) 



As in App. lA 31 we consider orbits with reflections from 
armchair and infinite mass edges to obtain for the sub- 
lattice part 



Tr (a a K s y a b K^ x ) 



Tr 



* a a? e ie ^ (1 + a aix >)° b (l + <wK 



(-l) N Ti 



(&a&b + Va,O aiX >O b O ai z 



= -4cos(6» a + 26 L )cos(6 b ) 

= (-l) NL 4cos(6 a )cos(6 b ) + 0(e), (A22) 
where we used Eqs. (|A16HA18|) and the geometric relation 
26 L = (N L + l)n + 0(e) mod 2tt . (A23) 



For the valley part we find that 



Tr 



t - 



= Tr [(KlY 

(-1) Nl+1 Tv(t ) for odd N ac , L 
Tr[exp(4iKA L T z )} for even A aCjL 

!)JV i+ i fo r odd iV aCii , 



= 2 

cos(4KAl) for even iV aCi £ . 
Thus for the whole pscudospin part, we have 



! cos(6> a ) cos(6 b ) 



(A24) 



(A25) 



(-1)^ cos(4A'A L ) for even N^ L , 

which can be generalized to include reflections from 
zigzag edges as in App. lA II yielding Eq. ([67]) . 



5. Trace y p p / for the universal conductance 
fluctuations 



As for the orbit pairs before, we consider quadruplets 
that hit only armchair and infinite mass type edges and 
consider the sublattice part first. Then we have 



Tr (a a K s y a a K*\ 



(A26) 



Tr 



o- o0 -fV^ CT * (1 + ow )a a (l + cr aiK )t 



N, 



= (-l) N yi N y +N -y 



x Tr 



= (-l) N yi N y+ N ^'A cos(^) 

x cos(6>; - 7 - 6 Y - [Ny + N y ]n/2) . 

We make use of the identities 

8^ -\- 6~,r = 26. 



(JV,,+ - JV,,_ - N r ,+ + 
w q ir + w r n + <D(si, e 2 ) 



N r ,-)l 



(A27) 
(A28) 



26> 7 = iV 7 7T + <p aiat ' a - (fl Xa a N mod 27T 

= (7V 7 + 1)tt + 6' a - 6 a mod 2tt (A29) 

and perform an analog calculation also for the second 
trace to obtain 



Tr \a a K$* a ief) Tr (a b K s p a b K s J 
16 cos(6> a ) cos(8' a ) cos(db) cos(6' b ) , 



(A30) 



if the orbit p (p 1 ) traverses the part q (r) in the same 
direction as 7 (7') [Fig.[8]a) and c)]. In the opposite case 
[Fig.|5]b) and d)], the result is 

Tr (a a K*a a Kty Tr (a b K s p a b K s J 
= (-l) NL 16cos(8 a )cos(6' a )cos(db)cos(8' b ) , (A31) 

with L = q + r^ 1 . Also for the valley part we have to 
distinguish the cases of parallel and antiparallel traversal 
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of the loop by the orbits p (p') and 7 (7'). We find that 
the valley part is zero for both cases, if N ac l is odd. For 
even N ac l we have for the parallel case 



Tr f K;k:}\ Tr (k;K v J) = Tr 2 (K v q K?) 



(A32) 



= 4cos 2 (2AW L + N im , L it /2) 
and for the antiparallel case 

Tr (k;k;J) Tr (k;K v J) = Tr {K* q K?) Tr (k^K^I 

= 4(-l) Ar "«. i cos 2 (2^A L + jV im , L tt/2) . (A33) 

Again this result can be generalized for zigzag edges 
yielding Eq. $SU\) . 



Appendix B: Loop density and spectral form factor 
in the presence of a magnetic flux 

1. Diagonal contribution to F (i) 

We consider a perpendicular magnetic field B, weak 
enough that wc can neglect effects on the pscudospin 
or the classical orbits, i. e. we only take into account 
Aharonov-Bohm like phases produced by the magnetic 
flux. An orbit 7 obtains a phase shift exp (i7r£M 7 /$o), 
where A 1 is the directed area enclosed by 7. For pairs of 
identical orbits, the phase shifts cancel, but for the pairs 
(7,7 _1 ) the resulting phase difference is 27r£?A Y /<£ > - For 
chaotic systems the distribution of the A 7 is known to be 
approximately Gaussian^ 



1 



^/2iraT~ 



exp 



2aT~, 



(Bl) 



with a system specific parameter a. Integrating the 
Aharonov-Bohm type phase shift over the enclosed areas 
with density Pa results in an exponential suppression of 
the contribution of the corresponding orbit pairs: 



00 

J dA 1 P A (A^,T 1 )e 2 ^ BA ^' ] ' a =e 



with 



Tr 



a 2 n n 

2n 2 a $ 2 ^ $ 2 ' 



(B2) 



(B3) 



Since the pairs of identical orbits are not affected, we 
obtain for the diagonal contribution to the spectral form 
factor 

F D (t) =t(l + e -2«VT^ ^ + e -tT H /T B ^ ^ (B4) 



which turns into Eq. (|30| for Tb — > 00, i. e. zero magnetic 
field. For very strong ac scattering, T ac — > 0, we re- 
cover the known formula for the GOE to GUE crossover 
Schrodinger billiard^. 



2. Loop contribution to F(t) 

For the loop contribution the procedure to incorporate 
a flux is similar. Here we have always a phase differ- 
ence exp (i2TrBA R /$o) between the orbits 7 and 7 X , so 
that the Tr integration in Eq. (|48|) is modified by an ex- 
ponential such as in Eq. (|B2j) . However, one has to be 
more careful here when taking into account orbit pairs 
like (7, 7* ), since they have the same propagation di- 
rection along R and the opposite direction along L in 
contrast to the pairs (7,7x)- Without magnetic field, 
this difference was irrelevant, but for finite B we have to 
distinguish these contributions. The T R integral has to 
be modified by a factor 



00 

J dA R P A [A R ,T R - T n 



is e 



2iriBA R /* 



-[T B -T min (e)]/T B 



(B5) 



for the pairs (7, 7x) anci by a factor 

00 

J dA L P A [A L ,T 7 - T R - T min (e)]e 2 ™ s ^/*° 



-[T 7 -T H -T min ( E )]/T B 



(B6) 



for the pairs (7,7x 1 )^ 2 -. Analog to our discussion in 
IIII B 41 we obtain for the loop contribution to the spectral 
form factor 



F L {t) 



-2t 2 



Tn 



-tT H /T B 



T b + Tac/2 



(B7) 



x 1 



— {e tTli 
T a( 



-2tT H /T ai 



For zero magnetic field, Tb — > oo, this formula turns into 
the former expression (|51[) . while for finite magnetic field 
and very strong ac scattering, T ac — > 0, wc obtain again 
the known expression for the GOE to GUE crossover in 
a usual Schrodinger billiard^: 



F L (t) « -2t 2 e~ tTH/TB 



(B8) 



Appendix C: Connection of Kubo formula and 
scattering matrix formalism: Fisher-Lee relation and 
Landauer formula 



First wc derive an equation that connects the Green's 
function of a scattering system to the elements of its scat- 
tering matrix via projection on the lead channels, i. e. a 
graphenc version of the Fisher-Lee relations^. To this 
end we define the wavefunction of a graphene lead in 
channel m as 



*m0«0 



(Cl) 



where x and y are the local coordinates of the lead such 
that the positive ir-axis points away from the scatter- 
ing system. The signs ± denote propagation in positive 
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(i. e. outgoing modes) and negative (i. e. incoming modes) 
x-direction. While the transverse wavefunctions ip m {y) 
of states with different m are orthogonal in the case of 
the Schrodinger equation without magnetic field, this is 
in general not the case for graphene. However, the cur- 
rent operator is still diagonal in this subspace (as for the 
Schrodinger equation with a magnetic field^l), namely 
for wavefunctions normalized to unit flux we have 



(C2) 



(C3) 



where C is the corresponding lead cross-section. To sec 
this we consider two eigenmodes m and m! that fulfill the 
Dirac equation 

{<y x k m + o y p y /h)\i/} m ) = k E \i' m ) , (C4) 

+ a y p y /h) = k E (ipm'\- (C5) 

We multiply the first equation with (ip' m \ from the left 
and the second with \ip m ) from the right and subtract 
the resulting equations to get 



(k m > - fe m )(Vwka#m) = 0. 



(C6) 



Then Eqs. ((C2J and Qg3\) follow directly: If 
k m i — k m 7^ 0, the matrix element has to vanish. Note 
that m m' does not necessarily imply k m / — k m =/= 0, 
e. g. for metallic armchair leads or if one considers 
spin degenerate modes. However in this case one can 
construct an orthogonal basis in the degenerate subspace 
so that the above still holds. 

Now consider a scattering system that is connected 
to an arbitrary number of leads. If an electron enters 
from lead /3 in the mode to', the wavefunction in the 
asymptotic region, i. e. far away from the scatterer, reads 



$(as) 



*m'(^)+Er!i-*+(^) ifxep, 



Tt {a(3 l^ + (x ) 



if x G a/(5, 
(C7) 

where the sums run over all propagating modes m in the 
corresponding lead. On the other hand the equation of 
motion for the retarded Green's function gives 

= J dV G(x, x')(k E - in ■ ^ x > + irj) ${x') , 

V =&(x-x') 

(C8) 

with V;,.' acting to the left on the second argument of 
G. Here V is a region in space that includes the scatterer 
completely and its boundary dV intersects the leads per- 
pendicularly in the asymptotic region. Integration by 
parts then leads to 



=i j dy' a G(x, x' a ) a Xa m , (x' a ) , 



(C9) 



where n x > is the unit vector normal to the boundary V 



at the point x' and a x 



We used that in the 



asymptotic region the retarded Green's function contains 
only outgoing waves and can be expanded as^i 



/ r OO — J_ j 



(CIO) 



Further we employed the current orthogonality relations 
(fC2|) and JC3|. Now we multiply both Eq. fC7]) and 
Eq. (|C9[) with the appropriate outgoing wavefunction in 
lead a and use again Eq. (|C2[) to project out the trans- 
mission amplitude 



(a/3) 



(Gil) 

= iv F J dy a J dy' p ^ m \x a )a Xa G(x a ,x' p) a X0 *~(a^) 

yielding the grap hene version of the Fisher-Lee relations 
(see also Ref.l73l). 

We can now derive the Landauer formula for graphene 
that expresses the conductance as a sum over transmis- 
sion probabilities of the individual propagating modes 
in the leads. To this end we note that in analogy to 
Eq. (jClpp . the advanced Green's function can be con- 
structed solely out of incoming modes in the asymptotic 
regional. Thus, for the retarded Green's function it fol- 
lows that 

G(x a ,x' p ) x ^°°J2^ a )dU4), (C12) 
in' 

and in combination with Eq. (jClpp we have 



G(x a ,x' ) Xq '^°° ^ f mm ^+,(x a )^J(x' p ) . (C13) 



Applying Eq. (|C2j) twice gives for the coefficients 

fmm' = ii>F tmm> ■ We insert this expansion into the linear 
response expression (1521) for the conductance and finally 
use Eq. (|C2[) again twice to obtain the Landauer formula 



5^ = xEic^i 2 = T Tr ( ttt )' 



(C14) 



where we defined t, 



In reverse, we can obtain Eq. (|52[) from the Landauer 
formula (|C14[) and the generalized Fisher-Lee relation 
(|CTTt using 



E / dy^ m ,{y)^,{y)v x ^ m {y) = ± 



(C15) 



which follows directly from Eq. (|C2[) . Inserting (|C11[) into 
the Landauer formula (|C14[) and using Eq. (IC15I) twice we 
get then Eqs. and JS5 
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In the same way, we further derive Eqs. (|88|) and (l89l) . Eq. (|C11|) . wc obtain 
that are relevant in the context of shot noise. Again with 

Tr(tthtf) = ^0 t* no t np t* mp =V 4 F dyi... dy s 

x ^(yi)<7x ,G(x 1 ,x 2 )<T Xf! ip~ (y 2 ) 

x ^t ] {y^) a x a G(x5,xe)<T X p ~>Pp (ye) 

x ipp Hy7) G\x s ,x 7 )(T Xa -0™(ys) • (C16) 

With the structure of the Green's function (|C13j) wc then 
use Eq. (|C15[) four times and obtain after renaming the 
integration variables expressions and 
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